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ABSTRACT 

In  this  paper  convergence  of  finite-dimensional  controllers  for  infinite-dimensional  sys¬ 
tems  designed  using  approximations  is  examined.  Stable  coprime  factorization  theory  is 
used  to  show  that  under  the  standard  assumptions  of  uniform  stabilizability/detectability, 
the  controllers  stabilize  the  original  system  for  large  enough  model  order.  The  controllers 
converge  uniformly  to  an  infinite-dimensional  controller,  as  does  the  closed  loop  response. 
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1  Introduction 

We  consider  semigroup  control  systems  on  a  Hilbert  space  X: 

x{i)  =  Ax(t)  +  Bu{t),  (1) 

y{t)  =  Cx{t), 

a:(0)  =  xq,  Xo  €  D{A)  C  X. 

The  operator  A  generates  a  strongly  continuous  semigroup  of  operators  ^(t)  on  X  so  A  is 
closed  with  domain  'D(A)  dense  in  X.  We  assume  bounded  control  and  observation,  and 
that  the  input  and  output  spaces  are  finite-dimensional:  B  G  B{R'^,X),C  G  B{X,Rp).  This 
control  system  will  often  be  abbreviated  (A,  B,  C).  For  control  functions  in  £2(0,  oo\  R^)  we 
interpret  the  solution  of  (1)  in  the  mild  sense: 

y{t)  =  CS{t)xo  + C  f  S{t  —  T)Bu{T)dT. 

Jo 

Further  details  can  be  found  in,  for  instance  [14]. 

The  equations  (1)  are  a  model  for  a  number  of  control  problems,  including  those  where 
the  system  dynamics  are  described  by  partial  diflferential  equations  and  hereditary  differential 
systems.  A  closed  form  solution  can  be  computed  only  in  the  simplest  of  situations.  In 
general,  it  is  necessary  to  use  an  numerical  approximation  to  the  semigroup  control  system 
(1)  in  order  to  simulate  the  response  of  the  system  .  This  approximation  will  typically  be  a 
system  of  n  ordinary  differential  equations,  which  we  write  as 

x{t)  =  Anx{t)  +  B„u{t),  (2) 

y{t)  =  CnX{t), 
a:(0)  =  aron- 

This  approximation  is  also  used  to  compute  controllers  for  the  original  system  (1).  Further 
details  on  the  approximation  scheme  will  be  presented  in  a  subsequent  section.  The  above 
finite-dimensional  system  will  often  be  abbreviated  as  (An,5n,C'„). 

In  this  paper,  we  are  concerned  with  finite-dimensional  controllers  of  the  form 

i(0  =  (A„  -  BM^t)  +  Fn{y{t)  -  (3) 

U„(t)  =  I<nZ{t). 

There  are  a  number  of  convergence  questions  associated  with  this  approach: 
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•  Does  the  controller  (3)  stabilize  (1)  for  large  enough  n? 

•  Do  the  controllers  (3)  converge  to  an  infinite-dimensional  controller? 

•  What  relation  does  the  closed  loop  performance  of  the  controller  (3)  with  Cn) 

have  to  its  implementation  with  the  original  system  {A,  B,  C)? 

Convergence  in  the  graph  topology  [18,  19]  of  the  approximating  systems  to  the  origi¬ 
nal  system  is  necessary  for  the  validity  of  an  approximation  scheme  as  a  basis  for  controller 
design.  Convergence  in  the  graph  topology  is  equivalent  to  convergence  of  a  sequence  of 
coprime  factors  (A^„,D„)  for  the  approximations  to  a  coprime  factorization  for  the  original 
system.  Without  such  convergence  it  cannot  be  concluded  that  a  controller  which  stabi¬ 
lizes  also  stabilizes  (y4n,5„,C„)  for  sufficiently  large  n,  or  that  the  closed  loop 

responses  converge.  In  fact,  in  the  absence  of  convergence  in  the  graph  topology,  the  set  of 
controllers  which  stabilize  both  the  infinite-dimensional  system  (A,B,C)  and  the  approxi¬ 
mations  {An,  Bn,Cn)  may  be  empty,  even  for  large  n. 

In  this  paper  this  approach  is  used  to  obtain  complete  answers  to  the  above  questions 
for  a  wide  class  of  systems.  A  complete  answer  to  the  problem  of  convergence  of  LQG  type 
controllers  is  given.  Although  the  graph  topology  is  the  fundamental  topology  underlying 
these  convergence  questions,  it  is  not  used  explicitly  in  this  paper.  All  results  are  derived 
directly  using  only  theory  on  Bezout,  or  stable  coprime  factors.  First,  some  background 
material  on  stable  coprime  factorizations  is  given. 

2  Stable  Coprime  Factorizations 

Suppose  a  system  of  the  form  (1)  with  xq  =  0  maps  inputs  in  €2(0,00,  to  outputs  in 
£2(0, 00,  R^),  and  that  furthermore,  there  is  a  maximum  ratio,  the  C2-gs\n  between  the  norm 
of  the  output  and  the  norm  of  the  input: 

II  2/llr:2(o,oo.Rp)  —  ”lll  “Ilc2(o,oo,/i”»)- 

Then  the  system  is  said  to  be  ^2-stable,  and  by  the  Paley- Weiner  Theorem,  the  Laplace 
transform  of  such  a  system  is  a  matrix  with  entries  in  Hoo  ■  Here  Ho©  indicates  the  Hardy 
space  of  functions  G(5)  which  are  analytic  in  the  right-half  plane  Re{s)  >  0  and  for  which 

sup  sup  lG(x  -f-  ju)\  <  00. 

>*'  x>0 

The  norm  of  a  function  in  Hoo  is 


II  <^'lloo  =  sup  sup  iG(x  -I- ju;)l. 

u;  T>0 
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We  denote  matrices  with  entries  in  "Hoo  by  M{  "Woo  )•  The  norm  of  a  function  in  M{  'Hex,  ) 
is  the  induced  matrix  norm 

II  ^11  oo  =  sup  sup  a  [G(a:  +  juj)] 

i>0  w 

where  a  denotes  the  largest  singular  value. 

The  Laplace  transform  of  a  linear  time-invariant  system  such  as  (1)  is  called  its 
transfer  function  and  is  given  by  CR{s\  A)B  where  R{s\  A)  indicates  the  resolvant  of  A. 
Suppose  G  is  the  transfer  function  of  a  given  system,  for  which  we  wish  to  design  a  controller 
with  transfer  function  H,  of  compatible  dimensions,  arranged  in  the  feedback  configuration 
shown  in  Figure  1. 


r2 


Figure  1.  Feedback  System 


The  2x2  transfer  matrix  A(G^H)  which  maps  the  pair  (ri,r2)  into  the  pair  (61,62)  is  given 
by 

'{I  +  GH)-'^  -G{I  +  HG)-^' 

A{G,H)  = 

H{I  +  GH)-^  {1  +  HG)-'^ 

The  feedback  system,  or  alternatively  the  pair  {G-,H),  is  said  to  be  externally  stable  if  each 
of  the  four  elements  in  the  above  matrix  belongs  to  the  set  S  of  stable  transfer  functions.  We 
could  define  stability  in  terms  of  the  transfer  matrix  from  (ri,r2)  to  (i/i,  j/2);  both  notions  of 
stability  are  equivalent  [19].  Definition  of  S  depends  upon  the  application.  Thus  the  closed 
loop  system  is  £2-stable  if  and  only  if  all  four  elements  belong  to  M(  Hex  )•  The  set  of  all 
plants  which  stabilize  G  is  written  S{G): 

S{G)  =  {H\A{G,H)eM{  Hex  )}• 

Note  that  the  present  definition  of  stability  is  symmetric  in  G  and  H.  Thus  G  stabilizes  H 
if  and  only  if  H  stabilizes  G. 

For  the  common  situation  where  the  system  is  already  stable  and  an  aim  of  controller 
design  is  to  improve  the  settling  time  of  the  system,  we  specify  a  real  number  a  >  0  which 
is  the  minimum  acceptable  stability  margin.  Then  a  system  is  said  to  be  input/output 
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<7  -stable  if  its  shifted  transfer  function  G{s  —  a)  is  in  M(  Hoo  )•  Equivalently,  we  replace 
"Hoo  by  the  algebra  of  functions  which  are  analytic  in  the  right  half  plane  Re{s)  >  —a, 
and  for  which 

sup  sup  \G{x  -f  jw)!  <  oo 

w  X>—(J 

with  corresponding  norm. 

Much  of  modern  control  theory  is  concerned  with  coprime  factorizations  of  systems. 
The  transfer  function  of  a  possibly  unstable  system  G  is  written  as  the  ratio  of  two  coprime 
stable  systems.  For  the  case  of  ^2-stability,  the  transfer  function  of  a  system  is  written  as 
G  =  ND~^  where  N,D  £  M{  Hoo  )  and  there  exists  X,Y  E  M{  Hoo  )  with 

X{s)N{s)  -1-  y(s)D(s)  =  I,  Ke{s)  >  0.  (4) 

{N,D)  is  called  a  right  coprime  factorization  (r.c.f.)  for  G.  Left  coprime  factorizations 
(l.c.f.’s)  are  defined  similarly.  {N,D)  is  a  l.c.f  for  G  if  G  =  D~^N  where  N,D  €  M{  Hoo  ) 
and  there  exists  X,V  €  M(  Hoo  )  with 

N{s)X{s)  -t-  Dis)Y{s)  =  /,  Re(s)  >  0.  (5) 

Every  system  which  is  described  by  a  system  of  linear  time-invariant  ordinary  differential 
equations  has  both  a  left-  and  a  right-coprime  factorization.  (This  is  a  consequence  of  the  fact 
that  the  transfer  functions  of  such  plants  are  composed  of  rational  functions.)  Furthermore, 
the  set  of  all  stabilizing  controllers  for  such  plants  may  be  described  in  terms  of  the  Youla 
parameterization]  a  controller  H  externally  stabilizes  G  if  and  only  if  it  can  be  written 

H  =  {Y  -  RN)-\X  +  RD),  \Y  -RN\i/^0,Re  M{Hoo)  (6) 

where  X^Y,N,D  are  as  defined  in  (4),  (5).  The  set  of  all  stabilizing  controllers  for  a  given 
system  G  are  parameterized  by  il ,  as  ranges  over  all  stable  systems.  In  other  words, 

S{G)  =  {{Y  -  RN)-\X  +  RD),  \Y  -  RN\^t),R£  M{Hoo)}-  (7) 

The  above  formulation  (7)  is  in  terms  of  left  coprime  factors  of  the  stabilizing  con¬ 
trollers.  The  same  family,  5(G)  may  also  be  written  in  terms  of  right  coprime  factors: 

5(G)  =  {{X  +  DQ){Y  -  NQ)-\  \Y  -NQ\^0,Qe  M{  Hoo  )}•  (8) 

Unfortunately,  more  general  systems,  which  do  not  have  rational  transfer  functions,  do  not 
necessarily  have  either  a  left-  or  a  right-coprime  factorization.  However,  any  system  with  a 
transfer  function  which  can  be  written  as  a  fraction  ND~^  D  G  M(  Hoo  ),  is  stabilizable 
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if  ajid  only  if  it  has  right  (and  left)  coprime  factorizations  [16]  .  For  systems  of  the  type  (1), 
the  theory  is  more  complete.  First,  some  definitions  are  required. 

Definition  1.1:  The  Co-semigroup  T{t)  is  stable  if  there  exist  constants  M  and 
a  >  0  such  that  |ir(f)||  <  Me~“‘  for  all  t  >  0. 

Definition  1.2  A  semigroup  control  system  (A,B,C)is  said  to  be  internally  stable  if 
the  semigroup  generated  by  A  is  stable  according  to  Definition  1.1  . 

Definition  1.3:  The  pair  (A,B)  is  stabilizable  if  there  exists  a  bounded  linear  oper¬ 
ator  K  :  X  —*  BJ^  such  that  A  —  BK  generates  a  stable  semigroup. 

Definition  1.4:  The  pair  (A,  C)  is  detectable  if  there  exists  a  bounded  linear  operator 
F  BP  —*  X  such  that  A  —  FC  generates  a  stable  semigroup. 

Definition  1.5:  The  system  (A,jB,C)  is  jointly  stabilizable/ detectable  if  {A,B)  is 
stabilizable  and  (A,  C)  is  detectable. 

The  extensions  to  a-stable,  <T-stabilizable  etc.  are  straightforward. 

Jointly  stabilizable/detectable  systems  (A,  5,C)  have  both  left  and  right  coprime 
factorizations  [5].  Let  K  be  such  that  A  —  BK  generates  a  stable  semigroup  and  F  such 
that  A  —  FC  generates  a  stable  semigroup.  Define 

N(s)  =  CB(s,A-  BK)B,  (9) 

D(s)  =  /  -  KB(s,  A  -  BK)B,  (10) 

A:(s)  =  KB{s,A-  FC)F, 

and 

Y{s)  =  1  +  I<B{s,  A  -  FC)B, 

XN  A  YD  =  1. 

{N,  D)  is  a  r.c.f  for  the  system  i.e.  CB{s-,  A)B  —  ND~^  [10]  and  the  Youla  parameterization 
describes  all  stabilizing  controllers  for  these  systems.  A  l.c.f.  can  be  defined  similarly,  using 
an  operator  F  such  that  A  —  FC  generates  a  stable  semigroup  .  Define 

ATs)  =  CBis,A~FC)B 

D{s)  =  I  -CB{s,  A-  FC)F. 

Then,  G(s)  -  D~^{s)N{s)  and 

NX  +  DY  =  I 
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for  some  X,Y  E  M{  Hoo  )• 

It  is  easy  to  show,  using  the  Hille-Yosida  Theorem,  that  internally  stable  control 
systems  (1)  are  externally  stable.  We  are  also  assured  that  external  stability  implies  internal 
exponential  stability: 

Theorem  2.1  [10]  A  jointly  stahilizable/detectable  semigroup  control  system  is  internally 
stable  if  and  only  if  it  is  externally  stable. 

This  equivalence  justifies  the  use  of  controller  design  techniques  based  on  system 
input/output  behaviour  for  infinite-dimensional  systems  of  the  form  (1).  The  discussion 
above  is  summarized  in  the  following  theorem. 

Theorem  2.2  Assume  (A,B,C)  is  jointly  stabilizable/ detectable  and  let  G  be  the  transfer 
function  of  the  system.  We  say  that  a  controller  with  transfer  function  H  stabilizes  G  if  the 
closed  loop  matrix  A{G,H)  €  M(  Tioo  )•  Then  H  stabilizes  the  system,  if  and  only  if ,  its 
transfer  function  can  be  written  using  the  Youla  parameterization  (7)  (or  (8  )  ). 

Also,  if  the  state-space  representations  of  the  systems  G  and  H  are  each  stabilizable 
and  detectable,  then  H  stabilizes  G,  according  to  the  definition  at  the  beginning  of  this 
section,  if  and  only  if  the  closed  loop  is  also  internally  stable. 

Theorem  2.3  [10]  Assume  {A,B,C)  is  jointly  stabilizable /detectable  and  that  H  is  a  con¬ 
troller  with  a  jointly  stabilizable /detectable  realization  {Ac,  Bc,Cc)-  The  closed  loop  system 
is  input/output  stable  if  and  only  if  it  is  internally  stable. 

The  following  is  by  now  well-known,  and  several  proofs  exist.  One  version,  which  uses 
coprime  factorizations,  and  which  first  appeared  in  [13]  is  given  here. 

Theorem  2.4  Every  jointly  stabilizable /detectable  system  (1)  is  stabilizable  by  a  finite¬ 
dimensional  controller. 

Proof:  Both  the  stable  factors  N,  D  have  well-defined  limits  at  infinity: 


lim  (t{N{s))  =  0, 


J?e(j)>0 


lim  a{D{s))  =  1. 

|«|-»oo 

Re(s)>0 


It  follows  (Mergelyan’s  Theorem)  that  each  factor  can  be  approximated  by  a  rational  element 
of  M(  'Hoo  ):  for  any  e  >  0  we  can  find  rational  Nn,  Dn  G  M(  Tioo  ),  with 

II  N{s)  -  Ar„(s)|L  <  e,  II  D{s)  -  T»„(s)|L  <  e. 
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For  sufficiently  small  e,  {Nn,Dn)  is  also  a  coprime  pair.  Let  (X,  K)  be  a  l.c.f  of  any 
finite-dimensional  controller  which  stabilizes  with  XNn  +  Y Dn  =  /•  Then,  if  e 

is  small  enough,  X  also  stabilizes  the  infinite-dimensional  system  since  X N  +  Y D  —  U 
where  U  has  an  inverse  in  M(  Tioo  )•  Rewriting,  we  have  X{NU~^)  -f-  Y{DU~^)  =  I.  Since 
{NU~^ ,DU~^)  is  a  r.c.f.  for  the  infinite-dimensional  system  (1)  ,  it  follows  that  every  jointly 
stabilizable/detectable  control  system  (1)  is  stabilizable  by  a  finite-dimensional  controller. 
□ 


3  Approximation  Scheme 

In  this  paper  we  are  interested  in  the  case  where  a  fmite-dimensional  controller  is  designed 
using  a  numerical  approximation  to  the  original  system. 

Suppose  we  have  a  sequence  of  finite-dimensional  subspaces  Xn  C  X  where  the  norm 
on  Xn  is  that  inherited  from  X.  Define  P„a;  as  the  orthogonal  projection  of  x  €.  X  onto  the 
finite-dimensional  subspace  Xn-  We  assume  the  following: 

(Al)  The  projection  operators  Pn  converge  strongly  to  the  identity  on  the  Hilbert  space  X. 
That  is,  for  all  x  €  X 

lim  11  P„x  —  x||  =  u. 

n-»oo 

For  each  the  approximating  system  is  (A„,J3„,C„): 

Bn  ■■=  Pr^B,  Cn  ■—  C|x„, 

and  An  is  an  approximation  to  A  which  satisfies  the  two  assumptions  (A2)  and  (A3)  below. 
Note  that  the  operators  An,Bn.Cn  and  the  semigroup  Tn{t)  generated  by  A„  are  operators 
on  Xn- 

A  core  C  of  a  closed  operator  A  is  a  linear  space  contained  in  the  domain  of  A  with 
the  property  that  the  set  of  elements  (x,Ax),  x  G  C  is  dense  in  the  graph  Q{A)  of  the 
operator  A  ([11],  pg.  166). 

(A2)  We  cissume  that  there  exists  a  core  C  for  A  such  that 

lirn  IjPnAx  —  A„P„x||  =  0,  for  ail  x  €  C.  (11) 

Such  an  approximation  scheme  is  said  to  be  consistent. 
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(A3)  We  will  further  assume  that  the  semigroups  T„(t)  generated  by  An  are  uniformly 
bounded,  that  is,  there  exist  real  numbers  N,M  >  1  and  k  such  that 

lir„(t)l|  <  Me*‘  for  all  n  >  N.  (12) 

(Uniform  boundedness  of  the  approximate  semigroups  is  generally  referred  to  as  “sta¬ 
bility”  in  the  numerical  analysis  literature.) 

Consistency  and  uniform  boundedness  are  sufficient  for  convergence  of  the  approxi¬ 
mation  [11,  14]  I'.e.,  for  all  £  >  0,<  >  0  and  for  all  a:  6  A  there  exists  N  such  that 

\[PnT{T)x  —  Tn(r)/’„x||  <  £  for  all  T  6  [0,t]  and  n  >  N.  (13) 

Assumptions  (Al)-(A3)  are  satisfied  by  typical  approximation  methods.  They  are 
sufficient  to  ensure  that  the  open  loop  response  of  the  systems  approximate  the  response  of 
G.  (Al)  is  of  course  redundant  in  that  it  is  implied  by  (13).  Assumption  (A2)  may  be  replaced 
by  a  requirement  of  strong  convergence  of  the  resolvants.  (Trotter-Kato  Approximation 
Theorem  eg.  [14]).  However,  these  assumptions  are  not  sufficient  to  ensure  convergence 
of  the  closed  loop  response.  An  example  of  a  scheme  which  approximates  the  open  loop 
response  of  a  system  but  is  not  satisfactory  as  a  basis  for  controller  design  is  given  in  [4]. 

The  appropriate  topology  in  which  to  establish  convergence  of  the  approximations  is 
the  graph  topology.  The  importance  of  the  graph  topology  in  controller  design  is  due  to  the 
following  result;  A  family  of  plants  G„  can  be  robustly  stabilized  by  a  compensator  H  which 
stabilizes  a  nominal  plant  G  if  and  only  if  G„  converges  to  G  in  the  graph  topology.  Fur¬ 
thermore,  in  this  case,  the  closed  loop  response  of  the  feedback  pair  A{Gn,H)  converges  to 
A(G,  H).  In  order  to  obtain  convergence  of  the  closed  loop  response  A(G„,  Hn)  to  A(G,  H) 
for  some  H  G  S{G),  Hn  must  also  converge  to  H  in  the  graph  topology.  The  results  given 
in  this  paper  can  be  shown  directly  using  the  coprime  factorization  theory  outlined  in  the 
previous  section.  The  graph  topology  is  not  directly  discussed  further  in  this  paper,  although 
it  is  the  topology  underlying  the  convergence  results.  The  interested  reader  is  referred  to  [19] 
for  an  introduction,  [20]  for  extensions  to  general  algebras,  [12]  for  its  application  to  approx¬ 
imation  of  semigroup  control  systems,  and  [6]  for  robust  stabilizablity  under  perturbations 
in  the  coprime  factors. 

Suppose  that  the  approximation  scheme  also  satisfies  assumptions  of  uniform  stabil- 
izability,  and  uniform  detectability: 

(A4)  If  the  original  system  is  stabilizable,  then  the  approximations  are  uniformly  stabi- 
lizable:  there  exists  a  uniformly  bounded  sequence  of  operators  {Kn}  such  that  for 
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sufficiently  large  the  semigroups  5„(t)  generated  by  An  —  Bnhn  are  uniformly 
bounded  by  Me““‘  for  some  M  >  0,a  >  0  and  all  n  >  N. 

(A5)  If  the  original  system  is  detectable,  then  the  approximations  are  uniformly  detectable; 
tnere  exists  a  uniformly  bounded  sequence  of  operators  {F,,}  such  that  for  sufficiently 
large  N,  the  semigroups  5„(t)  generated  by  An  —  F„Cn  are  uniformly  bounded  by 
for  some  A/  >  0,  Q  >  0  and  all  n  >  /V. 

Common  approximation  schemes  for  several  important  classes  of  systems  which  sat¬ 
isfy  these  assumptions  are  given  below.  Not  only  do  both  approximation  schemes  discussed 
below  satisfy  assumptions  (A1)-(A5),  there  exists  a  sequence  of  operators  I\„  wffiich  sat¬ 
isfy  assumption  (A4)  and  which  converge  strongly  to  an  operator  K  E  B{X,  U).  Similar 
convergence  exists  for  a  sequence  of  operators  F„  which  satisfy  assumption  (A5). 

3.1  Hereditary  Systems 

Consider  the  delay  functional  differential  equation 

x{t)  -  I  dT]{T)x{i  A  t)  +  Bu{t),  t>0  (14) 

J 

y{t)  =  Cx{t), 

i(0)  =  xq,  x{t)  =  —h<T<0 

where  x{t)  €  €  R’’,u{t)  €  F"*  and  B  and  C  are  matrices  of  appropriate  dimension. 

Also,  tj{t)  is  function  of  bounded  variation  taking  values  in  with  7/{r)  =  0  and  ?/(r) 

is  left  continuous  for  -h  <  t  <  0.  Defir.  ng  the  state-space  A'  =  F”  x  C2{  —  h,0\R^),  the 
hereditary  system  (14)  can  be  formulated  as  a  control  system  of  the  form  (1)  [2]. 

Define  the  finite-dimensional  subspaces  of  X  to  be 

Xn  =  |(<^1,02)  €  X;  M-r)  =  <  t  <  j  =  ....N  where  z,  G  . 

The  finite-dimensional  Galerkin  approximation  to  (14)  derived  using  these  subspaces  satisfies 
assumptions  (A1)-(A3)  [2,  15].  This  scheme  is  known  as  the  averaging  approximation  to  (14). 

In  [8]  it  is  shown  that  this  scheme  is  uniformly  stabilizable  and  detectable:  assump¬ 
tions  (A4)  and  (A5)  are  satisfied  .  It  is  also  shown  that  the  stronger  version  of  (A4)  and 
(A5)  holds.  For  a  bounded  operator  K  such  that  A  -  BK  generates  a  stable  semigroup,  we 
can  choose  =  A'P„,  A'„  satisfies  (A4)  and  lim„_oo  ||  A'„  -  A'||  =  0.  Similarly,  if  (14)  is 
detectable,  we  can  choose  a  sequence  F„  — »  F  where  A  —  FC  generates  a  stable  semigroup 
and  F„  is  a  sequence  which  satisfies  (A5). 
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3.2  Sectorial  Operators 

Define  A  through  a  continuous  sesquilinear  form  a  :  V  x  V  C 

{-A(f),ip)  ^  <f>  e  V{A),xl)  e  V.  (15) 

where  (•,•)  indicates  the  inner  product  on  X,  1/  is  a  Hilbert  space  ,  T>[A)  C  V  and  V  is 
densely  and  continuously  in  X.  In  order  to  avoid  confusion  with  the  norm  on  A^,  the  norm 
on  V  will  be  indicated  by  ||  Identify  A'’  with  its  dual  so  that  V’-^X  =  A'^V"'.  We 
assume  that  in  addition  to  (15),  a(,)  satisfies  Harding’s  inequality:  there  exists  k,c>  Q  such 
that 

a{u,u) k{u^u)  >  c||  (16) 

The  inequalities  (15)  and  (16)  guarantee  that  A  generates  an  analytic  semigroup  with  bound 

II  7^(011 

Further  details  may  be  found  in  [17].  Assume  that  the  approximating  subspaces  Xn  satisfy 
a  V'-approximation  property:  for  all  a:  €  there  exists  a  sequence  Xn  €  A^n  with 

(HI)  lim  II  Xn  -  xllv^  =  0. 

It  is  shown  in  [12]  that  this  class  of  semigroup  control  systems,  with  a  sequence  of  ap¬ 
proximating  subspaces  which  satisfy  assumption  (HI),  leads  to  a  sequence  of  approximating 
control  systems  which  satisfy  assumptions  (A1)-(A5)  of  the  previous  section.  Furthermore, 
as  in  section  (3.1)  we  can  choose  KniFn  so  that  /\„  K  and  —>■  F. 

Note  that  it  is  only  required  that  projections  onto  Xn  converge  in  the  V-norm.  A 
similar  result  in  [8]  requires  an  inverse  approximation  property: 

inf  II  /?(s;  A)z  -  x\\y  <  ei(n)||  2||  ^ 

x€  An 

inf  II  F{s;  A*  )2  -  x||v.  <  e2(n)||  ^11;^ 

where  »  0  as  A  — ►  oo.  These  conditions  are  stronger,  and  more  difficult 

to  verify  than  (HI).  Furthermore,  the  above  conditions  can  only  be  satisfied  if  R{s]  A)  is 
compact.  Problems  such  as  control  of  structural  vibrations  with  Kelvin- Voigt  damping  are 
thus  excluded. 
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4  Controller  Convergence 


It  will  now  be  shown  that  if  both  (A4)  and  (A5)  hold  in  addition  to  the  usual  (A1)-(A3), 
then  for  large  enough  model  order,  the  finite-dimensional  controllers  stabilize  the  semigroup 
control  system  (1).  If  a  stronger  version  of  (A4)  and  (A5)  holds,  then  the  controllers  con¬ 
verge  uniformly  (  in  M(  'Haa  )  to  an  infinite-dimensional  controller,  and  the  closed  loop 
performance  obtained  with  the  finite-dimensional  controllers  converges  to  the  closed  loop 
performance  obtained  with  the  infinite-dimensional  controller. 

Theorem  4.1  Assume  the  semigroup  control  system  (1)  is  jointly  stabilizable/detectable 
and  that  the  approximations  used  in  design  of  the  finite- dimensional  controllers  (3)  satisfy 
assumptions  (A1)-(A5).  Denote  the  transfer  functions  of  the  controllers  by  Hn  and  the 
transfer  function  of  (1)  by  G .  Assume,  in  addition,  that  the  feedback  operators  Kn  converge 
strongly  to  an  operator  K  G  B(X,U): 

lim  KnPnX  —  Kx,  X  ^  X. 

n— froo 

Then,  for  sufficiently  large  n,  Hn  stabilizes  the  infinite-dimensional  system  G.  For  such  n, 
the  operators 

Aon  ~ 

generate  stable  semigroups. 

Proof:  Define 

X„(s)  =  I<nR{s,  An  -  FnCn)Fn 

and 

y^(s)  —  I  KnR{s,  An  —  FnCn)Bn. 

It  is  well-known  (eg,  [19])  that  is  a  l.c.f.  for  the  finite-dimensional  controller  Hn., 

and  that 

XnNn  +  YuDn  =  I 

where  (A^„,D„)  is  the  r.c.f  of  defined  by  (9)  and  (10)  using  the  stabilizing  feedback 
operator  Kn¬ 
it  is  clear  from  the  assumptions  that 

II  ijv„  mil.  <  M 

for  some  constant  M.  Let  {N,D)  be  the  r.c.f.  of  G  defined  by  (9)  and  (10)  using  the  stabil¬ 
izing  feedback  controller  K.  It  was  shown  in  [12]  that  lim„_oo  =  V  and  lim„_oo  D„  =  D 


A  -BKn 

FnC  An  —  BnRn  ~  FnCn 
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and  so,  for  sufficiently  large  n,  say  n  >  Uq, 


It  follows  that 


II  [^n  Vn] 


N-Nn 

D-Dn 


<  1, 


n  >  Ho. 


=  XnN  +  YnD 


has  an  inverse  in  M{  Tioo  )  for  n  >  Uq. 

The  closed  loop  matrix  A{G,Hn)  can  be  written  in  terms  of  G  =  ND~^  and  Hn  = 

y;r'x„  as 

'l-NU-^Xn  -NU-^Yn' 

DU-^Xn  DU-^Yr,  ■ 

Since  A(G,  €  M{  Hoo  )  the  closed  loop  system  is,  by  definition  input/output  stable.  It 
is  trivial  to  show  that  Aon  is  the  infinitesimal  generator  for  the  closed  loop  system.  Internal 
stability  of  the  closed  loop  follows  from  the  fact  that  the  controller  is  stable  and  Theorem 
2.3.  Hence,  Aon  generates  an  exponentially  stable  semigroup.  Alternatively,  it  can  be  shown 
directly  that  the  closed  loop  system  is  stabilizable  and  detectable  and  Theorem  2.1  applied. 
□ 


The  dual  situation,  where  Kn  is  bounded  and  Fn  converges  strongly  to  an  operator 
F  ,  can  be  shown  identically,  using  convergence  of  a  sequence  of  l.c.f’s  of  the  systems  Gn 
and  boundedness  of  r.c.f.’s  of  the  compensators  //„•  We  also  have  the  following: 


Corollary  4.2  Let  a  be  as  in  assumptions  (A4)-(A5)  and  let  a  he  any  real  number, 
a  >  a  >0.  For  n  sufficiently  large,  the  operator  Aon  generates  an  a-stable  semigroup. 


Proof:  As  a  consequence  of  assumptions  (A4)  and  (A5)  we  can  replace  input/output  stability 
and  internal  stability  by  input/output  cr-stability  and  internal  cr-stability  in  the  previous 
theorem.  The  algebra  Tioo  is  replaced  by  the  shifted  algebra /fooa  • 

Theorem  4.3  Assume  that  the  approximation  scheme  satisfies  assumptions  (Al)-(A5)  and 
that  the  compensator  design  sequence  is  convergent: 


lim  KnPnX  =  Kx, 

n-^oo 


for  all  X  ^  X  and 


lim  F„y  =  Fy 

n—*oo 


for  all  y  £  IF  .  Let  H  denote  the  transfer  function  of  the  infinite- dimensional  controller 
whose  state-space  realization  is 


z{t)  =  {A-  BK)z{t)  +  F{u{t)  -  Cz{t)) 


(18) 
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y{t)  =  Kz{t). 

The  controller  transfer  functions  Hn  converge  in  M{  Tioo  )  i.e.  uniformly  in  s)  to  H. 


Proof:  We  first  show  that  l.c.f.  ’s  of  the  controller  converge: 

Jim  11  X  -  XnlL  =  II  ^  -  ^nlL  =  0- 

Define  Xn,Yn  as  in  Theorem  4.1  .  Denoting  the  semigroup  generated  by  An  -  F„C„  by 
Sfnit),  the  time-domain  representations  of  these  transfer  functions  are 

Xn{t)  :=  KnSFn{t)Fn, 

Yn{t)  :=  S{t)I  +  KnSFn{t)Bn 
where  S{t)  denotes  the  Dirac  delta  distribution.  Define 

Xit)  :=  KSF{t)F, 

Y{t)  :=  S{t)I  -h  KSF{t)B. 

The  proof  is  identical  to  that  used  to  show  convergence  of  coprime  factors  of  the  approxi¬ 
mating  systems  in  [12].  Due  to  convergence  of  Kn  to  K,  and  Fn  to  F,  we  can  show  that, 

lim  rilX„(O-X(t)lU„+pdf  =  0, 

n— »oo  Jq 

fOO 

lim  /  11  y;(0  -  =  o 

n-*oo  Jq 

where  jR'"+p  indicates  some  norm  on  operators  from  /?”*  to  B’’.  Convergence  in  this  norm 
implies  convergence  of  the  Laplace  transforms  of  these  functions  in  M(  7ioo  )•  Each  controller 
Bn  and  the  infinite-dimensional  controller  H  are  stable  (elements  of  M(  Hoo  ))•  This  implies 
that,  Yn,  Y  are  invertible  in  M(fH<x,)  and  since  lim„_oo  UK-  Y\\^  =  0, 

Jim  11  -  r-ML  =  0 

and  Hn  converges  to  H  in  norm.  (This  reflects  the  fact  that  the  uniform  topology  is  the 
restriction  of  the  graph  topology  to  stable  systems.)  Since  (X,  T)  is  a  l.c.f.  for  the  infinite¬ 
dimensional  controller  Lf,  the  result  follows.  □ 

Note  that  the  controllers  actually  converge  in  the  stronger  norm: 

lim  r  11  Hn{t)  -  =  0, 

n-*oo  Jq 

where  a  >  7  >  0  (a  as  in  (A4)  and  (A5))  although  this  is  not  used  here. 
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Corollary  4.4  The  closed  loop  operator  ^{Gn,  Hn)  converges  in  M{  TLoo  )  lo  A{G,H). 

We  can  write  the  closed  loop  operator  A{G,H)  in  terms  of  stable  coprime  factors  as 

and  similarly, 

The  result  now  follows  from  convergence  of  the  coprime  factors.  □ 

Thus  we  have  shown  uniform  convergence  of  the  controllers  to  an  infinite-dimensional 
controller,  and  similar  convergence  of  the  closed-loop  systems. 

4.1  LQG  Optimal  Control 

Commonly,  a  compensator  of  the  form  (18)  is  designed  in  two  steps,  through  construc¬ 
tion  of  an  optimal  regulator  and  an  optimal  detector.  Assume  that  (1)  is  jointly  stabiliz- 
able/detectable,  so  that  construction  of  an  internally  stable  closed  loop  is  possible. 

If  Q  is  a  self-adjoint  positive  semi-definite  operator  and  R  is  self-adjoint  positive 
definite,  then  the  optimal  solution  of 

/■oo 

J{u)=  /  {Qx{t),x{t))v  +  {Ru{t),u{t))dt 

Jo 

is  given  by  the  feedback  u{t)  =  —Kx{t)  where  K  :=  B*  and  11  is  the  unique,  non¬ 

negative,  self-adjoint  solution  of  the  algebraic  Riccati  equation, 

(  A*  n  -t-  HA  -  UBR-^  B*U  +  Q)z  =  0,  for  all  z  e  D{A).  (19) 

The  operator  A  —  BK  generates  a  stable  semigroup  on  X. 

Similarly,  we  may  obtain  the  observer  gain  F  by  solving  the  dual  problem 

AS  -1-  s  A*  -  SC*  r;'cs  +  =  0  (20) 

where  Qo  is  self-adjoint  positive  semi-definite  and  /?„  is  self-adjoint  positive  definite.  Setting 
F  =  E  C*  R~^  where  S  is  the  unique,  nonnegative,  self-adjoint  solution  to  the  above  Riccati 
equation,  A  —  FC  generates  a  stable  semigroup  on  X. 

Since  the  above  infinite-dimensional  Riccati  equations  cannot  be  solved  exactly,  ap¬ 
proximate  solutions  are  calculated  by  using  the  approximations  An,  Bn,  Cn-  Thus,  if  the  ap¬ 
proximations  are  stabilizable,  we  obtain  Kn  =  ^n  where  Hn  solves,  for  Qn  =  PnQPn, 

A*  n„  -f  IlnAn  —  YlnBnR  *  B*  Hn  +  Qn  =  0.  (21) 


I-NnXn  -NnYn 
DnXn  DnYn 


I -NX  -NY 
DX  DY 
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Similarly,  if  the  approximations  are  detectable,  the  operators  C*  Rg  ^  are  found  by 

solving,  for  Qon  =  PnQoPn, 

A„S„  +  <  -  Sn  C:  R:^Cn^n  +  Qon  =  0.  (22) 

The  theory  on  convergence  of  solutions  of  approximating  Riccati  equations  is  fairly 
complete,  and  we  can  state  the  following  which  is  a  consequence  of  [3]  and  the  theorems  in  the 
previous  section.  Convergence  of  the  solutions  requires  that  the  adjoint  semigroups  S*{t) 
converge  strongly,  uniformly  on  bounded  intervals,  to  S*{t).  This  assumption,  together  with 
(A1)-(A5)  implies  the  stronger  versions  of  (A4)  and  (A5)  required  for  Theorems  4. 1-4.3. 


Theorem  4.5  Assume  that  the  approximation  scheme  satisfies  the  usual  (A1)~(A3)  and  that 
it  is  uniformly  stabilizable  (A4)  and  uniformly  detectable  (A5).  Assume  also,  that  for  every 
X  E  X ,  the  adjoint  semigroups  S*  (t)  converge  strongly,  uniformly  on  bounded  intervals,  to 

5*  (t); 

lim  sup  II  SI  {t)P,,x  -  P„  S*  {t)x\\  =  0,  T  >  0.  (23) 

0<t<T 


Then: 


1.  The  controller  (3)  obtained  by  solving  the  finite-dimensional  ARE’s  (21)  and  (22)  to 
obtain  KniFn  stabilizes  the  infinite-dimensional  system  for  all  sufficiently  large  n. 

2.  This  closed  loop  system  is  internally  a -stable,  for  any  a  <  a. 

3.  The  controller  transfer  functions  converge  uniformly  (i.e.  in  M{  Hoc  ))  to  the  trans¬ 
fer  function  of  the  optimal  regulator/ observer  (18)  obtained  by  solving  the  infinite¬ 
dimensional  ARE’s. 


4.  The  closed  loop  operator  A{Gn,  Hn)  converges  uniformly  to  A{G,H). 


Proof:  It  has  been  shown  [3]  that  these  assumptions  are  sufficient  to  guarantee  that 
lim„_ooA'na:  =  Kx  and  that  \\mn-.oo  EnV  =  Fy  and  that  these  feedback  operators  yield 
uniformly  stable  semigroups  :  (A4)  and  (A5)  are  satisfied.  The  conclusions  now  follow 
immediately  from  Theorems  4. 1-4.3  .  □ 
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5  Conclusions 


A  similar  result  to  Theorem  3.1  and  Corollary  3.2  has  been  shown  in  [8].  However,  the  use  of 
coprime  factorizations  leads  to  a  considerably  shorter  proof.  Also,  the  results  in  [8]  depend 
on  several  eissumptions  that  are  not  required  here.  Unlike  [8],  the  spectrum  determined 
growth  assumption  not  required  at  any  point.  Furthermore,  both  classes  of  approximation 
schemes  discussed  in  [8]  have  the  property  that  the  resolvants  of  the  approximations  converge 
in  norm  to  the  original  resolvant: 

Jim  II  J!(A,/1„)-B(A,/1)||=0 

and  this  property  is  used  in  the  proofs.  This  restricts  the  class  of  systems  under  consider¬ 
ation  to  those  with  compact  resolvants,  and  the  approximation  scheme  to  those  with  this 
convergence  property.  As  demonstrated  here,  these  assumptions  are  not  necessary. 

The  results  here  showing  uniform  convergence  of  a  sequence  of  controllers,  and  of  the 
corresponding  closed-loop  systems,  have  not  been  previously  obtained.  Gibson  and  Adamian 
[7]  study  the  problem  of  convergence  of  LQG  controller  design  for  flexible  structures,  when 
the  controller  is  designed  using  approximations  to  the  solution  of  the  partial  differential 
equation  describing  the  vibrations.  Numerical  results  in  [7],  for  finite-element  approximations 
to  structural  vibrations  indicated  convergence  of  the  controller  transfer  functions,  but  the 
authors  were  unable  to  prove  Theorem  (4.3)  in  the  general  case. 

One  important  consequence  of  convergence  of  the  controller  transfer  functions,  is  that 
it  justifies  the  use  of  such  order  reduction  methods  as  Hankel  norm  reductions,  balanced  trun¬ 
cations,  and  coprime  factor  reduction.  The  order  of  the  approximation  (2)  could  be  increased 
until  the  controllers  (3)  have  converged  sufficiently  in  M(  Hoo  )•  The  numerical  results  in 
[7]  for  control  of  structural  vibrations  indicate  that  at  least  in  some  cases,  convergence  is 
obtained  fairly  quickly.  The  controller  order  can  then  be  reduced  using  a  standard  technique 
such  as  balanced  realizations. 
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